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This  note  lists  tables  computed  at  what  is  now  Numerical  Analysis 
Research „ University  of  California,  Lcs  Angeles  and  before  July  ls  195U 
was  the  Institute  for  Numerical  Analysis.  National  Bureau  cf  Standards.. 

Publication  of  tables  is  expensive  and  funds  for  the  publication 

! 

of  statistical  tables  are  extremely  hard  to  find.  It  is  therefore 
unlikely  that  many  of  the  tables  listed  telow  will  be  published.  Also, 
many  of  the  tables  are  at  present  not  ready  for  publication  and  some  of 
the  tables  are  too  special  to  warrant  publication.  It  seems  worthwhile 
to  report  the  existence  of  these  tables  and  indicate  their  present  state. 

The  following  list  contains  tables  connected  with  probability  and 
statistics  which  ha^e  been  computed  at  NAR  „ Some  of  these  tables  were 
computed  as  an  end  in  themselves , ethers  were  by-products  in  computations 
or  were  computed  because  the  machine  codes  were  available  from  other 
problems.  In  general,  the  tables  have  not  checked  beyond  the  accuracy 
required  for  the  purpose  for  which  the  table  was  intended.  It  seems 
worthwhile  to  list  tables  which  are  not  completely  verified  because 
they  may  simplify  the  checking  procedures  if  and  when  the  functions 
are  recomputed. 


■^The  preparation  of  this  table  war  sponsored  (in  part)  by  the 
Office  of  Naval  Research,  USN  and  the  School  of  Aviation  Medicine.  USAF. 
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The  number  of  dipits  given  in  the  tables  may  seem  excessive  to 
many  statisticians.  The  number  of  digits  used  is  a consequence  of 
the  fact  that  the  tables  were  computed  an  the  SWAC  (Bureau  of  Standards 
Western  Automatic  Computer).  This  machine  operates  on  a basic  number 
length  o.f  36  binary  digits  (10.8  decimal  digits)  and  there  is  no 
point  in  using  less  than  an  integral  multiple  of  the  basic  length  ands 
in  general,  the  tables  retain  most  of  the  digits  used  in  the  computation. 

Publication  appears  probable  for  the  following  three  tables  and 
they  are  therefore  not  listed  elsewhere  in  this  note. 

1.  11  Tables  of  the  bivariate  normal  distribution  function  and 

related  functions.51  Collated  by  the  National.  Bureau  of 
Standards.  The  introduction  is  by  G.  Blanch. 

The  functions  given  are 


2.  "Table  of  salvo  kill  probabilities  for  square  targets" 

Applied  Mathematics  Series  UU,  National  Bureau  of  Standards. 
Introduction  by  A,  D.  Hestenes. 

This  table  gives  values  of  the  function 
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3 o "Emrsir teal,  power  functions  for  non  parametric  two  sample 
tests  for  small  samples1'  D.  Teichroew.  Submitted  to  a 
technical  journal, 

Tl'is  paper  fives  the  emnirical  frequencies  of  all  possible  rankings 
which  are  obtained  when  a sample  of  m from  NfO-  i)  and  a sample  of  n 
tr cm  N(  ■ S , x)  are  ranked  in  order  of  size,  for  (m„  n)  = (3„  2)  (3,  3) 
(h,  2)  ana  Ur.  3^  and  various  values  of  S . 

Alx  the  tables  listed  below  exist  on  punched  cards.  In  addition, 
two  tabxes  mentioned  _n  section  II,  namely,  those  containing  y(p.»  A) 
and  Q (p),  have  been  multilithed  and  a limited  number  of  copies  are 
available  for  distribution. 

The  tables  are  listed  under  five  categories; 

I,  Tables  associated  with  the  normal  distribution. 

II.  Tables  associated  with  the  Gamma  distribution. 

HI.  Tables  associated  with  the  t-distribution. 

IV.  Tables  for  selecting  samples  from  certain  distributions. 

V.  Miscellaneous  tables. 
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I„  Tables  Associated  with  the"  Normal  Distribution 
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The  range  and  size  of  the  tables  is  as  follows! 


Decimal  PjLaces 


Function 

Range 

Tabulated 

Accurac 

f(x) 

x = -12 .00 ( .02)12 .00 

32 

* 2? 

[F(x)]k 

x ~ -12 .00( .02)12  .00 

k = 1(1)19 

32 

51  25 

H(xi  b) 

x - -12.00(.02)12.00 
b * 1(1)19 

32 

4 25 

^ ( a„  b) 

a5  b - 1(1)19 

32 

* 25 

E(x  . N) 

3 

j = 1(1) N 

N = 1(1)22 

21 

* 19 

E(x.2:,  N) 

j - 1(1)N 

N = 1(1)21 

21 

A 19 

E(x,  « x ,j  N) 
j-  J 

i,  J • 1(1)N 

N = 1(1)20 

21 

4 19 

k (Si  dk  ft) 

<?  - -3 .2(.1)0(  .01)6.1* 
d - 1(1)9 
ft  - 0(1)1* 

8 
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II,  Tables  Associated  with  the  Gamma  Distribution 
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Two  SWAC  routines  have  been  developed. 

1.  The  first  computes  p when  y is  given  for  X an  integer  by 
summing  the  series 


p(ys  A) 
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The  routine  is  also  used  to  compute  y for  a given  p by  inverse 
interpolation . 


V, 


2 0 The  second  routine  computes  y for  a given  p by  the  asymptotic 


series  (using  terms  up  to  and  including  Q^) 


y(p;  A)  = c2  + Q1(p)  c + Q2(p)  + 


Q,(p) 


where  A = c 


These  routines  have  been  used  to  compute  the  following  tabless 


Function 


p(y>  A) 


y(p>  A) 


2y(p>  A) 


Range 


Decimal  Places 
Tabulated 


Accuracy 


Q±(p) 


y -0(.5)  varying;  * 1*5.0  10 

A - 2(1)15 

p - .000( .001) .999  8 

A - 2(1)15,  20(10)50,  100 

p - .000( .001) .999  7 

2A  - 3(1)5(2)29 

p - .500(.00l).999  8 

i - 1(1)11 


III.  Tables  Associated  with  the  t distribution 


The  tables  are  concerned  with  solving  the  equation 
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for  t(pj  n)  when  p and  n are  given.  An  approximation  to  t(p;  n)  is 
obtained  by  summing  the  first  8 terms  of  the  asymptotic  series 
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The  tables  available  are  the  followings 


Function 


Range 


Decimal  Places 

Tabulated  Accuracy 


Vp) 


p - „500( .001) „999 
i * 1(1)8 


t(pj  n) 


p = 0500(o00l)999  8 

n -2(1)  16,  20,  25,  50,  100,  200 


IV  Tables  for  Selecting  Samples  from  Cei bain  Distributions 


One  method  of  selecting  samples  from  distributions  consists  of 


the  following  steps; 


1.  Compute,  9(k),  the  sum  of  k variates  uniformly  distributed 


on  (0,1) 


2o  Let  y - ZmQ  © 


The  ai  are  chosen  so  that  y has  the  required  distribution „ These 


coefficients  have  been  computed  for  the  case  where  k - 8 and  for  the 


following  distributions  ? 


1„  The  normal  distribution 


2 „ — | where  t has  a normal  distribution,  for  n = 50(1)200 


3„  The  Gamma  distribution,  for  A “ 2(1)  15,  20,  25,  50,  100, 


-i-  where  y has  a Gamma  distribution  for 
^ A - 2(1)  15,  20,  25,  50,  100, 


The  variates  — and  — are  used  to  generate  random  values  of 


\fn  y/j 


Wish art  matrices, 
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V.  Miscellaneous  Tables 

1„  Multinomial  coefficients „ 

This  table  gives  the  function 

(m+n+p+2)| 

”^~raTTupI 

for  all  combinations  of  m,  n,  and  p such  that  m + n + p 6 18. 

The  function  occurs  in  the  expression  for  the  covariances  of  order 
statistics . 

2.  Coefficients  for  curve  fitting  by  Chebyshev  polynomials. 

This  table  gives  10  decimal  digit  values  of 

cos  (2*.  + 1) 

for  n - 2(2)  12,  16,  18,  20,  25,  30,  40;  k - l(l)  n - 1,  &.  - 0(1)  n - 1. 
This  function  appears  in  the  curve  fitting  method  described,  for 
example,  in  "Tables  of  Chebyshev  Polynomials  Sn(x)  and  Cn(x)," 

National  Bureau  of  Standards,  Applied  Mathematics  Series  9. 

Introduction  by  Cornelius  Lanczos0 

3.  Tables  for  Probit  Analysis  with  Poisson  Error  Models. 

Three  functions  have  been  tabulated* 


) 


